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1 Introduction 



Thermodynamic Bethe Ansatz (TBA) is a method to calculate exactly the ground- 
state energy of an integrable quantum field theory in finite volume using its infinite 
volume scattering data [l] 1 The equations can be extended to excited states as well 
by analytical continuation [3,4]. 

The idea of the TBA is to exploit that the Euclidean partition function is domi- 
nated for large imaginary times by the groundstate energy. Calculating the partition 
function in the doubly Wick rotated (mirror) theory the imaginary time becomes the 
physical size which is taken to be large. Since the large volume spectrum is un- 
der control, the partition function can be evaluated in the saddle point approxima- 
tion which results in nonlinear integral equation for pseudo energies leading to an 
exact description of the ground state energy. Excited states on the complex (vol- 
ume/coupling) plane are connected to the groundstate which enables one to derive 
nonlinear integral equations for excited states as well. 

We start in Section 2 with a toy model containing one single particle with AdS 
dispersion relation and with scattering matrix which is not a function of the differ- 
ences of the momenta. Although this is a fictitious system it helps to introduce con- 
ceptual notions and steps needed to explain the TBA which is, in analogy, used in 
Section 3 to present the results for planar AdS/CFT. Finally, we give a guide to the 
literature in Section 4 and list some open problems. 

2 The concept of TBA: a toy model 

The application of the TBA method to solve completely the finite volume spectral 
problem is standard by now and follows the following steps. First the scattering 
theory has to be solved in infinite volume by determining the scattering matrix from 
its generic properties such as symmetry, unitarity, crossing relation. The poles of 
the scattering matrix lying in the physical strip are related to bound-states. These 
bound-states have to be mapped and their scattering matrices have to be determined 
from the constituents' scattering matrices. Then in the second step these scattering 
matrices can be used to describe the spectrum for large volume, which amounts to 
restrict the allowed particles' momenta via phase shifts and periodicity, and use the 
dispersion relation to express the energy in terms of the quantized momenta. This 
method sums up all power like corrections in the inverse of the volume and provides 
an asymptotical spectrum. The very same asymptotic description of the mirror theory 
is also needed as it can be used to calculate the exponentially small finite energy 
corrections from the partition function. Evaluating the Euclidean partition function 
for large imaginary times (large mirror volumes) in the saddle point approximation 

'The method has its origin in the work of Yang and Yang applied for spin chains and for the Bose 
gas with S interaction |2|. 
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provides integral equations describing the ground state energy exactly. Finally, these 
equations can be extended for excited states by analytical continuation. Now let us 
see how these steps are elaborated in the simplest setting. 

Infinite volume characteristics of the model 

We consider a toy model with one single particle type only. The dispersion relation 
is supposed to be the same as in the AdS/CFT correspondence 2 : 



The sine function indicates lattice behavior and restricts the momentum as 
p G [ — 7T, ?r] . The square root, however, has a relativistic origin. The theory is 
supposed to be integrable, thus multiparticle scattering matrices factorize into two 
particle scatterings. As relativistic invariance is not supposed the two particle S- 
matrix can depend separately on the two momenta S{p\,p2) and satisfies unitarity 
S{P\iP2)S{p2,P\) = 1 and crossing symmetry, which helps to fix it completely. We 
will not need its explicit form, but will suppose that in the p\ = pi = p particular 
case S(p,p) = — 1. 

Infinite volume characteristics of the mirror model 

The Euclidean version of the model is defined by analytically continuing in the time 
variable t = iy and considering space x and imaginary time y on an equal footing. 
The Euclidean theory so obtained can be considered as an analytical continuation of 
another theory, in which x serves as the analytically continued time x = H and y is 
the space coordinate. The theory defined in terms of y, z is called the mirror theory 
and its dispersion relation can be obtained by the same analytical continuation E = ip 
and p = iE which results in 



Contrary to the original theory the mirror model is not of the lattice type as its mo- 
mentum can take any real value As the scattering matrix is related via the 
reduction formula to the Euclidean correlator the mirror S-matrix is simply the ana- 
lytical continuation of the original scattering matrix: S(p\,p2)- 

Very large volume solution: asymptotic Bethe Ansatz for the model 

Let us put N particles in a large volume L subject to periodic boundary condition. 
Integrability ensures that the particle number is conserved and the particles' momenta 

2 The string tension is related to the 't Hooft coupling as 2ltg = VX 
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are not changed in the consecutive scatterings. The leading effect of the finite volume 
is the momentum quantization constraint: 

\=e^ L f\S(pj,p k ) (2.1) 

k:k+j 

which is called the Bethe Yang equation or asymptotic Bethe Ansatz (ABA) and fol- 
lows from the periodicity of the multiparticle wave function. Due to the sine function 
in the dispersion relation and the periodicity of p consistency of (2.1) requires L to 
take integer values only. 

Bound-states of the theory are manifested in the ABA as complex string-like so- 
lutions. Indeed, if the scattering matrix has a pole for Zm(p) > 0, then complex p 
solutions are also allowed in (2.1). If we take L very large with p\ rj ^ + iq then 
the rhs. of (2.1) would go to oo which should be compensated by another complex 
momentum, p2 ~ § — iq say, such that S{p\,p 2 ) exhibits a pole. The two particles 
with momenta p\ and p 2 form a bound-state with momentum p = p\+p%, energy 
E 2 {p) =E(pi)+E(p 2 ) and scattering matrix S 2 i(p,pj) = S(% + iq,pj)S(% -iq,pj). 
In general complex solutions built up from more particles are also allowed and they 
usually form a string-like pattern. Their dispersion relation and scattering matrices 
can be calculated by extending the method above, which is called the S-matrix boot- 
strap. 

Very large volume solution: ABA for the mirror model 

In the mirror model the considerations go along the same line as in the original theory. 
If we denote the mirror volume by R the ABA reads as 

l=e i * R HS(p j ,p k ) (2.2) 

k:k+j 

Since S(pi,p 2 ) lives in a different analytical domain than S(p\,p 2 ) its pole struc- 
ture can be also different. If it exhibits poles also at the proper location the mirror 
theory has also bound-states. Once bound-states exist we can calculate their disper- 
sion relation and scattering matrices from the bootstrap method. Suppose that the 
bound-states can be labeled with some charge Q, they have energy Eq{p) and their 
scattering matrix is SQjQ k (pj,pk)- The generic ABA valid for all the excitations (also 
for bound-states) will be 

\=e^ R Y\S Q] Q k {p b p k ) (2.3) 
Once these equations are solved the energy of the multiparticle state is 

7=1 

which describes the spectrum asymptotically for large volumes R. 
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(E,P) = (P y = iE,P x = P) 


(P y = P,P x = iE) = (E,P) 



Table 1: The relation between the original and the mirror model. 



Groundstate TBA equation from the partition function 

Let us come back to the original model and see how the exact groundstate energy 
can be determined in a finite volume L from the Euclidean partition function. We 
exploit the fact that the imaginary time evolution for large times, R, is dominated by 
the lowest energy state 

lim Z(L,R) = lim Tr( e - RH{L) ) = lim e - RE ^ L) + ... 

where the ellipsis represents terms exponentially suppressed in R. The same partition 
function can be determined alternatively, using the time evolution of the mirror theory 
which is generated by the mirror Hamiltonian H: 

Z{L,R) = Z(R,L) = Tv{e- Lfi ^) = 

n 

The relation between the original model and the mirror model is summarized in 
Table 1 . In switching to the mirror model we ensure that the volume goes to infinity 
(and not the imaginary time) where the spectrum is controlled by the ABA (2.3). 

In the large R limit the sum in the partition function is dominated by finite density 
particle states. Introducing the density of the particles (and bound-states) in momen- 
tum space (Pq(p) = |^) the energy can be expressed as 

E[p] =RZ [ d PPQ(P)E Q (p)= R L [ dup Q {u)E Q {u) 
Q J Q J 
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where for later convenience we reparametrized the momentum as p(u), momentum 
integrations go from — oo to oo. The quantization condition comes from taking the 
logarithm of the mirror ABA 

PA u i) + T, / du'{-ilogS QjQ/ ( U j,u'))p Q f(u') = 271^ (2.4) 
Q'J R 

where nj labels the quantized momentum pj whose charge is Qj. For a generic 
multiparticle state there are momenta pk which satisfy the same equation but which 
are not excited, not present in the system. They are called holes and their densities in 
the large volume limit is described by /5g. Clearly the densities of particles and holes 
are not independent they are connected by the thermodynamical limit of eq. (2.4) as 

d u p-27Z(p Q + p Q ) = -]T / duK Q g(u,u')pQi(u) =: -K QQ i-kp Q/ (2.5) 

& J 

where the kernel is defined as 

K QQ ,(u,u) = -id u logS Q& (u,u) 

The particle density itself does not characterize properly the states we sum over in the 
partition function. Indeed in a given interval (u,u + du) the occupied RpQ(u)du parti- 
cles can be distributed ( R Rp^fy" ") different ways leading to an entropy factor in 
the sum. Since in the large particle number limit the factorials can be approximated 
with the Stirling formula the partition function will take the form 

Z(L,R) = =£ f d [p Q ]e- LE ^\+s[PoM 

n Q J 

where the entropy factor is 

S[Pq,Pq]=R J du[(p Q + pQ)log(p Q + p Q )-p Q logp Q -p Q logp Q ] 

One can slightly generalize the partition function by adding a chemical potential term 
to the energy -LE Q [p Q ] -> \i Q [p Q ] - LE Q [p Q ] where ^Iq[Pq\ = RHq'Lq j dup Q (u). 
For fermions we take \Iq = in , while for bosons \Iq = 0. This extended partition 
function can be evaluated in the saddle point approximation. Taking into account the 
relation between 8pQ and 8pQ originating from the variation of (2.5) we obtain the 
minimizing equation in the so called pseudo energy £q = log ^ as 

du' 



/dii 
-^K Q , Q (u\u)\og(\+e 

=: -(log(l + e -^)*% Q )(") 
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Once the pseudo energies are determined the ground state energy in volume L can be 
obtained as 

E Q (L) = - E/ ^(duP) log(l + (2.6) 

The nonlinear integral equation which determines the pseudo energies is called the 
thermodynamic Bethe Ansatz (TBA) equation. Although it is not possible to solve 
it in general it provides an implicit exact description of the groundstate energy. This 
implicit solution is a starting point of a systematic large and small volume expansion 
and can be used to derive either functional relations for the pseudo energies or TBA 
equations for excited states by analytical continuation. 

Excited states by analytical continuation 

Here we start with bosonic theories without bound-states and suppose that by ana- 
lytically continuing in some parameter (say in the volume) we can reach all excited 
states. The way how excited states appear can be understood by analyzing the energy 
expression (2.6) integrated by parts 

E = J ^/5(«Klog(l+ e - e W) 

Let us suppose that in the analytical continuation singularities of type 1 +e~ e ("<) = 
appear. When we deform the contour their residue contributions give rise to 

£ = J ^^( M )log(l+e- e M) 

where we took into account the relation between the energy and the mirror momen- 
tum E(uj) = ip{uj). Taking the same analytical continuation in the equation for the 
pseudo energy we obtain 

e(u)=LE(u) + Y,logS(ui,u)- / ^K(w,u)log(\+e- e{w) ) 

. J 271 

Solving these equations iteratively for large L we can recognize that the 1 +e~ e ("<) = 
equations, which determine the positions of the singularities, coincide at leading 
order with the ABA equations (2.1). The subleading order calculation provides a 
universal formula for the leading finite size correction of multiparticle energy levels 
[5]. Alternatively for doing the analytical continuation one can think of the final 
result as choosing a different integration contour which surrounds the 1 + e~ e ("<) = 
singularities, and when we take the integration contour back to the real axis we pick 
up the above residue contributions. 
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Finally we note that if we have more species (labeled by Q) with diagonal scat- 
terings (like in the previous subsection) then a singularity in 1 +e~ e e, ("') = results 
in the equations 

E Q {u) = LE Q (u) + £log% e ( Mi -, u) - (log(l + e- £ Q' ) *K0 Q )(u) 

i 

whose solutions £q{u) and {«,•} have to be plugged into the energy formula 

£ = ££ G .("0-I / d ^d u p Q (u)\og{l+e-^) 

One has to be careful with such an analytical continuation in the presence of bound- 
states. Bound-states require pole singularities of the scattering matrices which usu- 
ally cross the integration contour in the analytical continuation and result in extra 
source terms. See the Lee- Yang model in the relativistic case [3,4] for example. 

3 TBA for planar AdS/CFT 

In this section we push forward the TBA program for planar AdS/CFT. The main 
difference compared to the previous discussion lies in the nondiagonal nature of the 
scattering matrix. There is a way, however, how we can profit from the previous 
diagonal results: the nondiagonal nature of any theory can be encoded into a diagonal 
theory but with auxiliary degrees of freedom. These auxiliary excitations do not 
contribute to the energy merely modifies the allowed momenta. Let us now follow 
the steps of Section 2. 

3.1 Infinite volume characteristics of the model 

The symmetry algebra of the theory has a factorized form: su{2\2) ®su{2\2). The 
fundamental particle called magnon transforms in the bifundamental representation 
whose S -matrix has the structure 

§n(pi,P2) = s {P^^Pi)Su{pi:P2)®Su{p\,P2) (3.1) 

where the matrix part S is fixed from its covariance under one copy of su(2\2) up 
to a scalar factor, which is determined from unitarity and crossing symmetry. The 
scattering matrix has simple poles corresponding to bound-states. There is an infinite 
tower of bound-states labeled by a positive integer charge Q. They transform under 
the tensor product of the atypical totally symmetric representations of the algebra and 
have dispersion relation 

Eq(p) = J Q 2 + 4g2 sin 2 | 
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3.2 Infinite volume characteristics of the mirror model 



As the mirror model is derived from the same Euclidean theory the fundamental parti- 
cles' scattering matrix is the analytical continuation of the scattering matrix (3.1 ). We 
are in a different analytical domain, however, and here different poles correspond to 
bound-states. These bound-states are also labeled by the charge Q but they transform 
under the atypical totally antisymmetric representations and have dispersion relation: 



E Q {p) = 2arcsinh Q g VP 2 + Q Z ) 



3.3 Very large volume solution: ABA for the model 

If we put N particles in a finite volume L the momenta of the particles will be quan- 
tized. The multiparticle wave function has to be periodic in each argument, that is 
when a particle transported along the cylinder it scatters on all other particles before 
arriving back to its initial position. In a diagonal theory this results in (2.1 ). In a non- 
diagonal theory, however, the multiparticle transfer matrix has to be diagonalized. 
This can be achieved by introducing new type of (magnonic) particles with vanish- 
ing dispersion relations and considering the original problem in terms of them as a 
diagonal scattering theory. 

Here we focus only on the charge Q = l sector of the theory. We have momentum 
carrying particles (»i) which scatter on each other as 
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1 



1 



S'uiPhPz) = S(p hP2 ) = ^ Tl-^°u. 



*J1 



where x ± (p) = ' 2 g ( 1 + yj 1 +4g 2 sin 2 |J and a represents the dressing phase. 

These particles are extended for each ,™(2|2) factor with two types of auxiliary parti- 
cles (y, oj), whose parameters are labeled by y G R and w G M. The auxiliary particles 
have trivial dispersion relations (their energy and momentum are zero) and scatter 
with the fundamental, momentum carrying ones as 



s?y(p>y) = -r^ y f- = ^pr l ; sSM = i 

Furthermore, they scatter on each other as 

Sn(wi,w 2 ) =S- 2 {w\ -w 2 ) ; Sy\(y,w) = Si(v(y) -w) ; 5^(yi,y 2 ) = l 



3 The index 1 in »i refers to the charge of the particle. This particle is a first member of an infinite 
series of bound-states labeled by »q. Similarly we will meet particles of type o N and \>m- 
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where v(y) = y +y 1 and we introduced a useful function S n (v — w) = V _ W _L ■ Any 

scattering matrix can be extended by unitarity to the opposite order of their particle 
types/arguments: S(i,j)S(j,i) = 1. 

In formulating the ABA equations for the full theory we have to take into account 
the two su(2\2) factors and that they commute. The ABA equation for the momentum 
carrying particles reads as 

i=e^ L Y\s\\{p hPk ) n n^yipjtf) 

k:k+j a=l,2/=l 

where N' is the number of fundamental and N y a the number of y type particles, while 
the a = 1 , 2 index refers to the two su(2\2) factors. Since the two factors commute the 
ABA equations for the auxiliary particles with rapidities y 1 ' 2 and w 1 - 2 can be written 

as 

ft w,Pk)n%<yj,wf) = i = n w>>?) n 5 n»n 

k:k+j 1=1 k:k+j Vl+k 

Not all solutions of the ABA equations correspond to single trace operators as the 
level matching/zero momentum condition has to be fulfilled £ y pj = 0. The theory 
contains also bound-states which can be determined from the singularity structure 
of the scattering matrices. Since from the TBA point of view only the bound-state 
spectrum of the mirror theory is relevant we will focus only on them. 

3.4 Very large volume solution: ABA for the mirror model 

In the case of the mirror theory the fundamental scattering matrix is the analytical 
continuation of the original one p — >■ p. As a result the ABA will be the analytical 
continuation, too 

1 = e^ R J] SZiPM EI TltfyiPjtf) ( 3 - 2 ) 

k:k+j a=l,2/=l 
k:k+j 1=1 

1 = n^K'^)ft 5 nK'<) (3-4) 

t.k^j Ujtk 

There are some differences compared to the original ABA. First the domain of p G R 
is different compared to p G [—71, n] and the total mirror momentum does not need to 
vanish. Then, as we are in the mirror theory, the way how is expressed in terms 
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of p is also different: x 




(p-0 



1 + T 1 J ■ Additionally, in the calculations 



of the ground state energy the sectors with antiperiodic fermions are relevant and this 
is manifested in a minus sign in the middle equation. The possible bound-states and 
their ABA equations are the subject of the next section. Let us note that usually in the 
literature instead of (3.4) its inverse is considered as this will lead to positive particle 
densities in the thermodynamic limit. 

3.5 Exact groundstate energy: TBA 

In this section we derive TBA integral equations for the groundstate energy in finite 
volume R. We treat the theory as if it were diagonal with the scattering matrices 
specified above. First we analyze whether this "diagonal" theory has bound-states by 
analyzing the thermodynamic behavior of the equations and calculate the scattering 
matrices of the bound-states, the so called strings. They are special complex solu- 
tions of the ABA equations and they all contribute to the partition function which 
determines the ground state energy. Then we use the canonical procedure to derive 
coupled integral equations for the pseudo energies in a raw form, finally, using identi- 
ties between the scattering matrices originating from the symmetry, we rewrite them 
in a simplified form and analyze simple excited states. 

3.5.1 String hypothesis for the mirror model 

The string hypothesis is similar to closing the S-matrix bootstrap program, that is 
to identify all particles (including bound-states) of the theory and to determine their 
scattering matrices. Let us premise that we will find bound-states of three infinite 
types (»q, >m, °n) for Q,M,N G N, and also of a finite type y§ particle with 5 € {±}. 
They can be arranged in the two dimensional lattice shown in Figure 1 . Let us see 
how they arise from the ABA equations. 

In the following we put R and all particle numbers large (keeping their ratio fi- 
nite) and analyze the ABA one by one. Let us first note the reality properties of the 
equations. Unitarity of the mirror scattering matrix implies that the y roots come in 
complex conjugated pairs y,- = (yj 1 )* or lie on the unit circle y = (y -1 )*, similarly 
the roots w come in complex conjugated pairs w, = w*or are real. 

•q particles 

In looking for momentum bound-states we rewrite the scattering matrix in (3.2) as 



S'niPuPi) 
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where the rapidity is introduced as =i ± + ^. As R is very large complex 

values for u\ with positive imaginary part are allowed. In this case the lhs. of (3.2) 
for 7 = 1 diverges so there should be another u say U2 that goes to Mi — j. If ui 
still has a positive imaginary part then by the same argument there should be another 
// say W3 which goes to u 2 — j. Applying this procedure we arrive at a string of Q 

roots K + (G-l)f,w+(fi-3)|,...,w-(e-3)| ! M-(G-l)f or shortly u Q+l _ 2 j = 
u+ i(Q + 1 — 2j) I where j = 1 , . . . , Q. (Clearly the Q=l string is the original particle 
itself.) The scattering of the Q-string with any other particle of type (.), label i and 
rapidity q is 

2 

SQi(u,q) = n 5 *! : ( M 2+l-2M) = S' Q {q,u)- 1 

7=1 

Although naively the scattering matrices seem to depend on the parameters and 
such a way the bound-state scattering matrix depends on its constituents, this is not 
the case when we take into account the contributions of the dressing phase as was 
shown in [6]. 

The auxiliary particles exist for both su(2\T) factors. Here we focus only on one 
of them and omit to write out its index. 

yg particles 

Let us analyze (3.3). If we suppose that the number of momentum carrying particles 
N' goes to infinity then 

n; r if |y 7 -|<l 

n S$(yj,Pk) -> { ±1 if M = i (3-5) 

<:<-'./ [ 00 if \yj\ > 1 

In the middle case y roots lying on the unit circle are allowed. As the scattering 
matrix 5^,° (y,w) has a difference form in the variable v(y) = y + y _1 we might use 
the parameter v instead of y. The inverse of the relation, however, is not unique. 
Defining y_(v) = i(v — iy/4 — v 2 ) with the branch cuts running from ±00 to ±2 we 
can describe any y with 3m(y) < for v G [—2, 2]. Clearly y+(v) = y_ (v) _1 describes 
the other 3m(y) > case and in the scattering matrices 5^* which depends on y, and 
not on v, we have to specify which root is taken. As a consequence we have two 
types of y particles yg with 5 = ± and the scattering matrices split as S y y \ (y, q) — > 
S y s \(ys{v),q) =:S y s \(v,q). 

>m particles 

If |yi| < 1 in (3.5) then the rhs. of (3.3) goes to zero which has to be compensated 
by a w\ root which goes to vi — 1 = y\ +y\ l — L - But then taking the ABA for w\ 
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means that the rhs. of (3.4) will diverge which has to be compensated by a root 
V2 = w\ — |. If the corresponding y2 satisfies |j2| > 1 then (3.3) is consistent with 
(3.5) and reality requires y\ = {y^ )*, w\ = w\. The three roots yi -H- vi = v+| 
and wi = v and v — ^ = V2 j2 form an M = 1 string which we denote by >i. In 
the case when |j2| < 1 then we have to repeat the same arguments for y2 leading 
to W2 and yj, and so on. Finally we arrive at the notion of a >m string. It consists 
of 2M y particles with yj = (yZj)* and M o particles with synchronized parameters 

wm+i-2; = v + (M + 1 - 2 j) l - and y } -»• v sign(7 - )(M+2 _2j) = v + sign(y') (M + 2 - 2/) | 
for j = 1, . .. ,M. The composite scattering matrix of the >m particle with all other 
particles is simply the product of the scatterings of its each individual constituents 

M+l M M-l 

= EI S -i( v M+2-2j,q)Y[ s xi( w M+i-2j,q) II ti+i( v u-2j,q)=S£[(i,v)~ 1 

7=1 7 7=1 

o w particles 

Suppose we have a large number of y particles and that w\ has a positive imagi- 
nary part. Then the first factor of the rhs. of (3.4) will go to zero which has to be 
compensated by a root W2 = W\ — — . If 3m(w2) < then we obtain a o 2 string. 
In the opposite case we repeat to previous argumentation leading to an /V string 
Wjv+i-2; = w + {N + 1 — 2j)^. Clearly a single w is just a oj string. The scatter- 
ing of the N string with any other particle is 

N 

S°m i w ' = II S wi ( w n+ 1-2 j, i) 

7=1 

Scattering matrices 

Summarizing, the mirror AdS theory in the thermodynamic limit could be replaced 
by a diagonal theory having constituents of infinite type (•,>, o) and index Q,M,N 
for Q,M,N G N, and also of finite type y particles with 8 6 {±}. See also Figure 1. 
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Table 2: Scattering matrices of the various particles 
For the readers convenience we summarize the scattering matrices in Table 2 
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The scattering matrices are unitary SySji = 1 and their explicit forms are 
Sq' q ,(u,u') = SQQi(u-u')L Q& (u,u'y 2 

Q-l 

Sqq'(u-u') = SQ + Q'(u-u')SQi-Q(u-u')Y[SQi-Q + 2j(li-u') 2 

;=1 

Q Q' 1 



Zqq(u,u') = Y\Y\o(u Q+ i- 2 j,UQ + i-2k) — 
j=lk=\ 1 



x{u Q -2j)x(u Q , +1 _ lk ) 

1 

x{l<Q + 2-2jH"Q'-2k) 



where uj = u + j~ and we reparametrized the momentum carrying particles in terms 
of the rapidity via the function x(u) = \{u — iy/4 — u 2 ). Recall also that S n (u — w) = 



"" - . The other matrix elements are 

U- n — w 



c .> , \ X(U- Q )-X(V M ) x(U-q)-x(V- M ) x{uq) n 

Jqm{ u , v ) = — ( \ 7 n 7 \ 7 V-—, C I I ?>M-Q-2j{U,V) 

* x(uq)-x(v M ) x(uq)-x{V- M ) X(U-q) f = \ 



x(u-q) — x(v) s j x(u 



x(uq)—x(v) S y x{U-q) 

S^ M i(u,u) = Smm'(u — u) = Sm M ,(u,u) 1 
5 ms("' v ) = s m{u-v) = S°J s (u,v) 
The ABA equations then have a generic form 

k a=l,2 I m n 

where . can be any type of •,>, o,y but only the • particles have nonvanishing energy 
Eq and momentum Pq{u) = g(x{u-o) —x(uq)) + iQ. The parameter F denotes the 
fermion number. We also indicated the contributions of the two su(2\2) factors. The 
energy of such a multiparticle state having Nq of Q k particles is 

E{px,...,p k )=Y i E Qk (p k ) 

k 

Let us note that the ABA equations for the auxiliary particles can be inverted with- 
out changing their physical meaning. Taking the inverse of (3.4) is equivalent to 
redefining simultaneously the scattering matrices S°^ s — > (S° y s )~ l and Sjj — > (Sji) -1 . 
Actually these are the equations used in the literature as they give positive particle 
densities in the thermodynamic limit. 



3.5.2 Raw TBA equations 

Suppose now that we would like to describe the groundstate energy in the AdS system 
in volume L. In doing so we follow the steps presented in Section 2 to evaluate the 
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partition function for large mirror sizes. We introduce densities of particles (strings) 
Pq{u), Pm{u), Pn{u) for u G R and p y s (u) for u G [—2,2] and the analogous densities 
of holes p — > p . They are restricted via the logarithm of the ABA which contains the 
logarithmic derivatives of the scattering matrices 

Kjj,(u,u') = -id u logSjj,(u,u') 

Clearly Kjj,(u,u') / —Kj,j(u',u) as the scattering matrices are not of the difference 
type. (Keeping in mind how we obtained the string solutions the densities are natu- 
rally ordered Pg S> p y 3> P^Pm-) Then we introduce the entropy factors for the den- 
sities, in chemical potential for fermions and calculate the saddle point of the func- 
tional integral. This results in integral equations for the pseudo energies £g, e^, e^, e y s 
as follows 

s' Q = LE Q - log(l +e- £ Q>) *K* Q : Q - log(l +e- e ") *K& Q - log(l +e~^) *K y ' Q 

where in the contributions of the >m and y§ we have to sum for the contributions of 
the two su(2\2) factors (which we omitted to write out). The remaining equations are 
valid separately for the two su(2\2) factors separately: 

e» M = -log(l+e- e J)^-log(l+e- e «')*^M-log(l+^)^M 
e° N = log(l+e- e ~0*^w + log(l+^)*^; 

e y s = -log(l+,- e a)^-log(l+e- e -)^-log(l+e- e -)^ + /7r 
Once these equations are solved the groundstate energy can be obtained as 

E o( L )=-if ^ G iog(i+,- e e) 

Finally we note that we replaced the magnonic ABA for the particle type o N with 
its inverse and made the corresponding change in the scattering matrices to ensure the 
positivity of the magnonic densities p^. It effectively changed the sign of the related 
kernels. 



3.5.3 Simplified TBA equations and Y-system 

In this subsection using identities among the TBA kernels we bring the equations in 
to a universal local form. This means that pseudo energies can be drawn in a two 
dimensional lattice, such that only neighboring sites couple to each other with the 
following universal kernel 

sl MN = S M n ~{K+l)mf > s ( u ) = 4co f h g^ (3 - 6) 
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where Imn = §m+i,n + §m-\,n and (K+l) M l N *(K NL + 8nl) = 8ml- To simplify the 
notation let us introduce the following Y functions 

Y- Q =e- £ h ■ Y M = e- £ - ; Y£ = e e ° N ; Y y § = e 5 ^ 

Clearly we have two copies for Y^j a ,Y^' a ,Yl ,a . (To conform with the literature we 
inverted the ABA equations for >m and y_). Acting with the operator (3.6) on these 
inverted TBA equations and using kernel identities like (K + 1) mn *Kn = s8m,i we 
arrive at their simplified, universal form 

1 1 + Y y 

logy^ = iog(i+y^ +1 )*5-/ MM /iog(i + — )*s+<5 M ,iiog- — 

i + yZ 

logF^ = I NN ilog(l+Y^)-ks + 8 Nj i\og- r *s 

1 + T^v 

where in the convolution * we integrate over the interval [—2,2] only. The other 
equations do not behave so nicely. 

log^ = -/ GG ,log(l + -^ r )* 5 + log(l+^ 1 )*^ + log(l+F^ 2 1 )* 5 ; <2>1 

io g y; = -iog(i + -^)* 5 +(io g (i+yZ' 1 )(i+rZ' 2 ))*^-A*s 

where A vanishes on the interval [—2,2] whose explicit form can be found in [7 ]. The 
equation for the y particles are simpler in the original form 

8 log Y y s = — log( 1 + Yq) *Kq S + log l -± Y f *K M + in 

These equations for Yg are not in a local form. However, acting with the inverse 
of s they can be brought into such form. The operator s~ l acts as (f * s~ l )(u) = 
f(u + i — iO) +f(u — | + iO) and involves the analytical continuation of the functions. 
It has a large null space, thus when acting on the equation information is lost: 

logFZ**- 1 = log(l +Y{) +log(l +yf) -log(l + -L) 

The advantage of defining s^ 1 in the above manner is that it uses the analytically 
continued values of the Y functions on the rapidity torus only. If we continue them 
across the cuts by using (f*s~ l )(u) = f(u + 1 — iO) +f(u — | — iO) = f + (u) +f~ (u) 
then the term A disappears, but the Y functions have to be extended to an infinite 
genus Riemann surface. On this surface the Y-system has the universal form 



(1 +yjy,M+i)(l + yjy,M-i) 
(1 +^V_lJif)(l +^v+i,m) 



I N,M I N,M~ 7i , v-1 \7l i v-1 ^ ' 
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Figure 1: Y- system for planar AdS/CFT. F_ is denoted by G while Y + by ©. 

where the 7Y,M indices live on a two dimensional integral lattice. In our situation 
the identification can be drawn on Figure 1, which explicitly reads as Yq = Fg,o, 
Yli a = Y M+UVa , Y°' a = Y 1MN+1) , Y y _> a = Fi, Va and Y* a = F 2 , Vq2 where Vi = 1 and 
v 2 = -1. 

3.5.4 Excited states by analytical continuation 

Here we focus on the TBA equations for excited states in the sh sector for small 
coupling. This sector contains particles of type »i only and have ABA: 

\=e^ L Y\S\l{p u , Pl ) 

These equations are asymptotic only and the exact system of TBA equations is re- 
quired to describe the energy of the multiparticle state exactly. As the vacuum is 
a BPS state it has vanishing energy and its analytical continuation cannot describe 
excited states. Alternatively we choose an integration contour, such that when it is 
taken back to the real axis the residue of a singularity of the form 1 + e~ e '( pk ^ = is 
picked up resulting in additional source terms in the raw equations as: 

£g log stg ; £M^L lo gS*MO;>") ; lo § s 's(Pj> u ) 

j j j 

Once the new system of TBA equations are solved the pseudo energies Sq have to be 
plugged into the energy formula: 

OO /> T 

= / *"d u p Q \og(l + e- e 'a) 

k Q=\J A71 

to obtain the energy of the multiparticle system. 

We can rewrite the TBA equations in terms of the F functions into their simplified 
form. They satisfy the same F-system relations (3.7 ) but with a different asymptotical 
behavior. There is a systematical asymptotical expansion of the F-system, which 
reproduces both the ABA and the leading Liischer correction of these multiparticle 
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states. This is valid for weak coupling g — > (or large sizes) and it is very nontrivial 
to follow the analytical behavior of the Y functions as one increases the coupling. The 
ABA solution itself suggests, that additional 1+7 = singularities could appear and 
then the TBA equations have to be modified by additional source terms. These source 
terms ensure the analytical behavior of the energy around these singular points. 

4 Guide to the literature 

Here we list the representative papers where the various parts of the TBA program 
were developed. 

The idea that the TBA program can be applied in the planar AdS/CFT setting was 
presented in [9]. The infinite volume scattering description of theory can be found in 
chapters [ 10 11]. The ABA equations for the planar AdS/CFT model was conjectured 
in 1 12] (and thoroughly discussed in chapters [13. 14]), while the analogous ABA 
for the mirror model was described in [15]. As the color structure (su(2\2)) of the 
scattering matrix is the same as that of the Hubbard model, the Hubbard TBA solution 
can be adopted [16]. This results in the string hypothesis which was formulated 
explicitly in [17]. The standard procedure leads to raw TBA equations, which were 
developed in [ 19, 18, 20]. The simplified form of the TBA equations was presented 
in [7] and the Y-system relations, presented previously in [21 ], were derived in [ 19, 
18,20]. In doing this the analytical properties of the dressing phase [6,20,22] had to 
be investigated. In the AdS/CFT context the volume of the integrable system has to 
be an integer, which can be seen also on the groundstate TBA [23]. 

Although we obtained the Y-system from the ground-state TBA equations, in 
principle, it follows from the hidden PSU(2,2\4) symmetry of the model. An in- 
dependent alternative approach based on this symmetry is the subject of the next 
Chapter in this volume [24]. 

The Y-system plays a crucial role in describing excited states. As it is related 
to the symmetry of the model [24-26] it is the same for each state. What makes 
the difference is the asymptotical and analytical behavior of the Y-functions. The 
analytical properties of the Y-functions was thoroughly analyzed in [19, 28, 23, 8 ]. 
Based on the solution of the Y-system of the 0(4) model [27] the authors of [21] 
identified the large volume solution in terms of the transfer matrices of the ABA [ 14 ]. 
This helps to derive excited states TBA equations for the sl2 sector, which was done 
in [20,28]. The excited state TBA equations provide an exact description of the 
given state and they were used in the Konishi case, [29, 30], to analyze numerically 
the behavior of the energy for large coupling. The results are summarized in Figure 
2 4 , see also [31 ]. It was further shown in [28] how to modify these excited state TBA 
equations if a 1 + Y = singularity appears in the analytical continuation in g. 

4 We thank the authors of 1 3 1 1 for borrowing their figure. 
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Konishi state 




Figure 2: Numerical solution of the excited TBA equations for the Konishi 
state [29 , 30 1 . 



The weak coupling limit of the Y-system equations can be compared to the ABA 
[14] and Liischer type correction [32]. The leading order behavior is built in the 
asymptotic solution [21] of the Y-f unction, but the next to leading one provides a 
stringent test of the excited states TBA equations, which was performed numerically 
for the Konishi operator in [ 33] and analytically at next to leading order in [ 34] . Later 
this analytical calculation was extended to describe the next to leading order Liischer 
correction of generic twist two states [35] in [36]. 

The strong coupling limit of the F-system for a finite density of string particles 
was analyzed in [37], where a complete agreement with the one loop string energies 
including all exponential finite size corrections has been found. The functional Y- 
system equations were encoded into simpler Q functions in [38,25, 31 1. 



v 

£ a,s 
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AF 


BFT 


GKKV 


%)(") 




Yq(u) 


Y Q (u) 


Y. Q (u') 


Ym+\,\{u) 




y mU u ) 


Y V: m(u) 








Y N \v(it) 


Y Wi n(u) 


Yo N+ M') 


%(") 


Yl(u) 


-Y+(u), 


Yy(u) 


Ye(u>) 


Yi,i(u) 


Yl(u) 


-YI l (u), 


Y r (u) 


Y®(u') 



Table 3: Relating the Y-functions to those in the literature, where u' = gu. 



Let us mention, how our TBA equations are related to those in the literature. We 
summarized the relation between the various conventions for half of the Y-system 
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in Table 2 as the other half is trivially related, see also [8]. Under this replacement 
our simplified equations are equivalent to AF [7], while the raw equations to BFT 
[18], except for the chemical potentials of [ 1 8 ] . In comparing to GKKV [20] the 
indentification is not enough. Comparing our kernel Kqq, to the one K, Q , Ql in [20] we 
observe a slight difference. This is irrelevant, however, for excited states satisfying 
the level matching/zero momentum condition 5 

The AdSs/CFT^ correspondence has a brother theory, the AcISa/CFTt, duality 
[39], where the TBA program has been developed in an analogous way. The ABA 
together with the string hypothesis of the mirror theory lead to ground state TBA 
equations and Y-system relations in [40,41] and extend the previously conjectured 
Y-system proposal of [21 ]. This program is further elaborated in [41] by additionally 
determining excited states TBA equations and comparing them to the asymptotic 
solution of the Y functions [21] and to the quasi classical string spectrum. 

Finally, let us list some open problems. 

There are two disagreeing string theory calculations ( [42] and [43]) for the 
anomalous dimension of the Konishi state. Additionally, the numerical solution of 
the TBA equations for large couplings [29, 30] provides a third result, and calls for 
improvements both the string theory and the TBA sides. On the string theory side 
it could be a pure spinor calculation, while on the TBA side one should analyze the 
analytical behavior of the Y-system and check whether, with increasing g, a singu- 
larity of type 1 +Y = indeed appears, as the asymptotic solution suggests [28]. In 
principle the effect of such singularities is to make the coupling dependence of the 
energies analytical, but it has to be established concretely. 

The anomalous dimensions of twist operators in the planar limit can be described 
by integral equations derived directly from the ABA [44]. It would be nice to see, 
how the exact excited TBA equations reduce to these equations in the large spin limit. 

The analytical comparision of the excited state TBA equations to the next to 
leading order Liischer corrections [34, 36] tested explicitly only the > part of the Y- 
system. A next to next to leading order analysis could test the o part as well. 

The excited states TBA equations are coupled nonlinear integral equations for 
infinite unknowns. An ideal system of equations should contain finite unknowns only, 
and could be developed in analogy to [27,45] by exploiting the result of [38,25,31 ]. 
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Note added in proof 

After this review chapter was finished three string theory calculations based on differ- 
ent methods determined the strong coupling expansion of the anomalous dimension 
of the Konishi operator [46 -48 ]. All agreed with each other and with the strong cou- 
pling expansion of the TBA equation [29, 30]. This gives a strong support not only 
for the correctness of the TBA equations but also for the integrability approach to 
planar AdS/CFT 
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